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ARATIONS 
POWHLU_R3 = H_floating ** unsigned integer power routine 
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+ 
FACILITY: Language support Library, Mathematics division 
ABSTRACT: 


set? module contains OTSSPOWHLU_R3, a procedure which raises an 
H_floating value to an unsigned” integer power. 
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ENVIRONMENT: Runs at any access mode, AST Reentrant 
AUTHOR: Steven B. Lionel, CREATION DATE: 22-JUL-1981 
MODIFIED BY: 

1-001 - Original. SBL 22-JUL-1981 
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¢$ -SBTTL DECLARATIONS 

48 : LIBRARY MACRO CALLS: 

0° SSSDEF : System error codes 

1 SCHFDEF ; Define condition handler symbols. 

¢ SSFDEF ; Define stack frame symbols. 

SPSLDEF ; Define program status longword 

: 3 symbols. 

: : EXTERNAL DECLARATIONS: 

8° -DSABL GBL ; Force all external symbols to be declared 
59 ~EXTRN MTHSK_UNDEXP ; Undefined exponentiation 

60 ~-EXTRN MTHSK_FLOOVEMAT ; Floating overflow in Math Library 
61 -EXTRN MTHSK FLOUNDMAT ; Floating underflow in Math Library 
6 -EXTRN MTHSSSIGNAL ; Math error routine 
64 : MACROS: 
65 ; 
? 3 NONE 
68 : EQUATED SYMBOLS: 
20 : NONE 
, : OWN STORAGE: 
ri : NONE 
0 : PSECT DECLARATIONS: 
78° .PSECT _OTSSCODE PIC, USR, CON, REL, LCL, SHR, - 
79 EXE, RD, NOWRT, LONG 
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-SBTTL OTSSPOWHLU_R3 = H_floating ** unsigned integer power routine | 


+ 
; FUNCTIONAL DESCRIPTION: 


OTSSPOWHLU - H floating ** 
1884 OTSSPOWMLU. R3 - 


Raise an Phare dey value to an uns |aned integer power giving 
a H_floating result in registers RO-R3. 


The result is given by: 


base exponent result 

any >0 product (base * 2**i) where i is each 
Non=zero bit position in exponent 

>0 = 0 1.0 

= 0 = 8 Undefined exponentation 

< 0 = 1.0 
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; CALLING SEQUENCE: 
result.wh.v = OTSSPOWHLU (base.rh.v, exponent.rlu.v) 
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; FORMAL PARAMETERS: 


00000004 base = 4 ; H floating base. Note that this is passed by 
immediate value in four argument List positions, 
which is a violation of the calling standard. 


This is allowed for language support procedures. 


00000014 exponent = 20 ; Unsigned longword integer exponent 


IMPLICIT INPUTS: 
NONE 
IMPLICIT OUTPUTS: 
NONE 

ROUTINE VALUE: 
The base raised to the exponent's power. This is returned in 
registers RO-R3, which is a violation of the calling standard, 
but is permitted for language support procedures. 

SIDE EFFECTS: 

May signal: 
MTH$_FLOOVEMAT - floating overflow in Math Library 
MTHS_FLOUNDMAT - floating underflow in Math Library, only 
if the caller has FU enabled. 
n 


MTHS$_UNDEXP - undefined exponentiation, if both base and 
exponent are zero. 
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8 | 
O1FC 13) «ENTRY OTSSPOWHLU_R3, “M<R2,R3,R4,R5,R6,R7,R8> ; Entry point 
9 128 VAB ate HANDLER, (FP) 3; Translate exceptions to MTH$ errors 
B 6 144 BISPSW gPSLsA_FU ; Enable floating underflow detection 
70FD A 145 VH 3 RO-R3 = initial result 
70FD 148 MOVH base tAP), RG : aa R = base 
08 1 14 MOVL exponent (AP) R8 : exponent 
1 1 198 BNEQ $ i f not zero 
ratR 19 14 TSTH kr 3 ase also zero? 
1 oF 131 BEQL UNDEF INED ; 06s error. Otherwise, answer is 
04 OO1EF 1 : RET : With result of 1.0 
a. 
OiF 155: Do the first iteration here so that we can clear the high-order exponent 
001" 136 : bit. Afterwards, that bit will be zero so we don't have to worry about 
O1F 158 i- 
Sol 180 BEGIN 
00 =€5 borF 161 BBCC #0, R8, 10$ ; Skip initial multiply 
70FD 00 ; 19¢ MOVH RO ; New partial result 
F 9C 8 : 197 10$: ROTL #-1, R8, RB : she 7 low exponent bit and clear 
13 60 . es BEQL DONE ; ise that att? If so, we're done. 
aes 
OO2E 168 Y Each time we get here we know that there's at least one more exponent bit 
33 : 199 R left, so square our current power of the base. 
OOF TN op 
64FD $6 138 "  MULH2 R4, R4 3; Square base. This may overflow or 
00 174 ; underflow. In either case, the 
003 175 : final result would also overflow 
003 176 3 or underflow. Our exception 
bos 177 ; handler will catch these cases. 
E9 00 178 BLBC R8, NEXT_BIT 3; See if this exponent bit is on. 
0 179 : If so, multiply base by product. 
64FD 0035 180 MULH2 R4, RO ; This too may overflow or underflow, 
039 #8 181 3 which will be caught by the handler. 
ae § NEXT_BIT: 
3 >. ASHL R8, RB ; Get next bit of exponent. 
1 ; : : BNEQ LooP 3; Keep going if not all done. 
oe 1 § DONE : : 
04 4 1 RET ; Return with answer in RO-R3. 
1 iB. 
rt} 139 + Undefined exponentation error - 0«*0 
41 192° 
41 198 UNDEF INED: 
9A 41 194 Roveet #ATHSK une XP, -(SP) ; jodigate undefined exponentiation. 
79 0045 195 ASHQ #15, #T, : reserved floating 
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SIGNAL_ERROR 
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MOVZBL #MTHSK_FLOOVEMAT, -(SP) ; 
ASHQ #15, #T, RO ; 


operand 


s@ 
; EXC_HANDLER unwinds here if we got a floating overflow. 


pieet ing overflow in Math Library 
RO-R3 = reserved floating 
operand 


+ 

EXC_HANDLER unwinds here if we got a floating underflow. 

If our caller enabled FU, then signal FLOUNDMAT, with a result of zero. 
Otherwise quietly return zero. 


LRH 
BBC 


RO : 
#PSLSV_FU, SFSW_SAVE_PSW(FP), DONE 


SIGNAL ~ERROR 


Initial result 
Just return if 


: caller disabled underflow. 
MOVZBL #MTHSK_FLOUNDMAT, -(SP) ; Floating underflow in Math Library 


3% 

; Signal a MTH$ error whose code has already been pushed on the stack. 

3; RO-R3 has either a reserved operand or a zero, depending on the 
exception type. 


wow Sete 


IGNAL_ERROR: 
CLRL 
CALLS 


RET 


(FP) 
#1, G*MTHS$SIGNAL 


Cancel our condition handler 
convert to 32-bit condition code 
and signal a MTH$ error. t 
second argument to show that this 
is a CALL entry point. 

peturn to caller with result in 
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6F 0 ie per forming the exponent iat on. 
; } ; then cone nue execution at OVERFLOW, as Be STH _FLOOVEMAT 
6F rk : Otherwise resignal. 
gt Seg 
GF 246 EXC_HANDLER: 
0000 Bot re «WORD “M<> 
50 oS AC D O71 i MOVL CHF SL eat rte RO 
AO D 75 50 TSTL CHFSL_MCH_DEPTH(RO) 
16 is 078 51 BNEQ RESIGRA a 
50 04 AC D O7A 26 MOVL CHFSL_SIGARGLST(AP), RO 
51 04 AO dO 07 5 MOVL FS$L_SIG_NAME(RO), R1 
04B4 8F 35 a BB 23 Hay Goo, OvERrL PLTOVE 
04C4 BF 51 B1 089 2$ CMPW Rist att FE TUND. F 
ii 13 OBE 5 BEQL foto _UNDERFLOW ~ 
389 58 RESIGNAL: 
50 0918 8F 3c }©60090 59 MOVZWL #SSS$_RESIGNAL, RO 
O* b098 508 ts 
0096 6 GOTO_OVERFLOW: 
51 60 00 8298 6 MOVL CHF$L_SIG ARGS ( RO), R1 
FC AO41 «BI AF ODE 0099 64 MOVAL B*OVERFLOQ, -4(RO)ER1) 
09 11 44 66 BRB CONTINUE 
O0A1 68 GOTO_UNDERFLOW: 
51 60 DO OQOA1 69 MOVL CHF$L_SIG_ARGS(RO), R1 
FC A041 B2 AF DE epee 4 MOVAL B*UNDERF LOW, -4(ROSCR1] 
oe are 3 BRE CONTINUE 
OOAA 74 CONTINUE: 
50 01 DO OOAA 75 MOVL #SS$_CONTINUE, RO 
04 OOAD 6 
OAE 7 
OOAE 78 
OOAE 79 END 
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Ext “geese = This condition pengter 9 gets. seperot if we got an exception 


overflow, 
to be 


loatin 


Similarly, continue at UNDERFLOW if an-underflow was seen. 


Entry point 


Get mechanism arguments list 
At depth zero? 

If not, resignal 

Get signal arguments List 
Get signal name 

Overflow fault? 

If so, continue at OVERFLOW 
Underflow fault? 

If so, cent tnue at UNDERFLOW 


Resignal exception 
Return to VMS exception dispatcher 


; Get number of signal arguments 


Move address of OVERFLOW routine 
to PC in signal argument List 
Continue execution 


Get number of signal arguments 
Move address of UNDERFLOW routine 
to PC in signal argument List 
Continue execution 


Continue at PC in etenet List 
Return to VMS exception dispatcher 


End of module OTSS$POWHLU 
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3.43 = 04 | 
HPs 1F R 02 
at L_MCHARGLST = | 
CHFSL_MCH_DEPTH = 
CHF $L-SIGARGLST = 4 
CHFSL_SIG_AR = 0 
CHFSL-S1G_NAME = 04 
CONTINUE AA R 8 
DONE 40 R 
EXC_HANDLER OOOOO6F R 0 
EXPONENT = 000 4 
GOTO_OVERFLOW 96 R 8 
GOTO-UNDERF LOW Al R 
LOOP™ 0000002E R 0 
MTHS$$S 1GNAL aeneeneer =X =O 
MTHS$K_FLOOVEMAT aeeeeeee =X 00 
THSK—F LOUNDMA eeeneeee X 00 
MTHSK_UNDEXP eeneeeee =X =O 
EXT BIT 00000039 f 
OTSSPOWHLU_R3 00000000 RG 0 
VERF LOW 00000040 02 
PSLSM_FU = si $44 0 
PSL$V_FU = 00000006 
RESIGNAL 00000090 02 
SF$W_SAVE_PSW = 00000004 
SIGNAL_ERROR 00000065 02 
85$. CONTINUE = 444 1 
SS$_FLTOVE _F = 00000484 
SS$_FLTUNDF = 000004C4 
SS$-RESIGNAL = 00000918 
UNDEF INED 00000041 R 02 
UNDERF LOW 00000059 R 02 
ee + 
! Psect synopsis ! 
tone nw eee meee eee + 
PSECT name Allocation PSECT No. Attributes 
: $ 00000000 < 0.) 00 ¢ QO.) NOPIC USR CON ABS LCL wosie NOEXE NORD wount wove BYTE 
$ 00000000 ( 0. 83 ( 1.) NOPIC USR CON ABS LCL NOSHR EXE’ RD OVEC BYTE 
OTSSCODE OOOOOOAE (¢ 174.) 02 PIC USR CON REL LCL SHR’ EXE RD NOURT NOVEE LONG 
} ; Performance indicators ! 
Phase Page faults CPU Time Elapsed Time 
Initialization 35 0: le 0:00:01 .98 
Command processing 115 0:00:00.5 : B98. 0 
Pass 206 0:00: 9) 00:13. 
syebo} table sort 0 :00: 9. 7 :00:01. 
Pas 61 :01.1 :00:02. 
sae E table output :00. :00:00.11 
Psect synopsis output 700. :00:00.05 
Cross-reference output 700. 0:00:00.00 
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| 
Assembler run totals 427 00:00:07.70 00:00:24.70 | 
Ihe york ing (83 Limit was 1200 pages. 
1 artes, (52 pages) of virtual memory were used to buffer the intermediate 
a re were 30 pages of symbol peste, space dicing. 18-0 d : hold 508 non-local ane Tiel symbols. 
¢/ source AS were read in Pass 1, produ bject records in Pass 
pages of virtual memory were used to deti ne 940 mechs. 
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Macro Library name Macros defined 
“S25SSDUAZB:(SYSLIBISTARLET.MLB;200 SS ; 
562 GETS were required to define 7 macros. 

There were no errors, warnings or information messages. 


MACRO/ENABLE=SUPPRESSION/DI SABLE=(GLOBAL , TRACEBACK) /LIS=LI1S$:OTSPOWHLU/OBJ=OBJ$:OTSPOWHLU MSRC$:OTSPOWHLU/UPDATE=(ENH$:OTSPOWHLU) 
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